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We calculate the scattering amplitudes of Goldstone bosons off the pseudoscalar D-mesons in 
unitarized heavy meson chiral approach. The low energy constants appearing in 0(p 2 ) chiral La- 
grangian are determined by fitting lattice simulations on S-wave scattering lengths. Z>s (2317) is 
obtained as a bound state in (S, I) = (1,0) DK channel. Possible bound states or resonance states 
in other channels are investigated as well. The quark mass dependence of the mass and binding 
energy of D* (2317) is also investigated, which indicates predominately DK molecular nature. 
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■ I. INTRODUCTION 

t> ■ 

CN| ' In the last decade, the discovery of many narrow resonances with open charm open a new chapter in hadronic 
spectroscopy. Especially, the £>* (2317) discovered by the BaBar Collaboration [T| and D s i(2460) by the CLEO 
Collaboration have inspired heated discussions both experimentally and theoretically. Moreover, the Belle Col- 
O ^ laboration recently reported a broad + charmed meson with mass and width being m D *o — 2308 ± 60MeV and 

T D *o — 276 ± 99MeV, respectively Meanwhile, the FOCUS Collaboration reported a broad + charmed meson 
. with mass and width being m D *o = 2407 ± 56MeV and r D .o = 240 ± 114MeV, respectively Although consistent 
fIZJ ' with each other within the errors, it is still in dispute whether they are the same particle 

Possible interpretations of -D* (2317) include normal cs state [6j, four-quark state Q, hadron molecular state 8|, 
' etc. To distinguish composite from elementary particles, different methods were proposed such as pole counting 9j, 
scattering length and effective range [l3| • As emphasized in a series of paper [TJ E2] , quark mass dependence of a 
C*~) ' state can also provide important information on its nature. 

Effective field theories(EFTs) have been proven very successful in studying low energy hadron physics In light 
meson sector, chiral perturbation theory is an expansion in powers of external momenta and masses of Goldstone 
bosons [3 El- In high energy region or for large quark masses, chiral amplitudes violate unitarity severely. In 
addition, chiral expansion up to a given finite order does not contain resonance or bound state, which may modify 
the results of physical variables from perturbation theory significantly. Therefore, unitarized model was introduced 
to high energy region, in which lower lying scalar and vector resonances can be dynamically generated. Following the 
same spirit, heavy meson chiral perturbation theory (HMChPT) was proposed [16H18| . and unitarization method was 
applied to some phcnomcnological analysis |19h23| . 

Lattice gauge theory is another powerful tool to study strong interactions. Lattice simulations are usually performed 
at unphysical quark masses, or equivalently at larger pion masses. Recently, lattice results for the charmed meson-light 
hadron scattering lengths are given at several chosen values of M-^/F^ [24j. These progresses can be used to make 
up the lack of experimental data on scattering processes. These lattice data can be used to determine the low energy 
constants in perturbative scattering amplitudes [25j |. 

D* (2317) as well as other possible charmed particles were investigated with the unitarized heavy meson chiral 
approach by studying the scattering lengths of charmed mesons and Goldstone bosons in Ref . [H, HH . The quark 
mass dependence of the poles has an interesting behavior and provides a good way to understand the structure of the 
obtained poles. In their calculation, the large Nc approximation and the mass and width of Z?* (2317) as input are 
used to determine the low energy constants. In this paper, we use the similar approach to reinvestigate the charmed 
mesons scattering off light mesons. The difference between our treatment and theirs is that we do not apply the large 
Nc approximation. This is because in the real world, Nc is 3. Moreover, in the large Nc limit, the mass and width 
of the particle can be quite different from the real particle [26|-[28j . Therefore, three additional low energy constants 
(LECs) ho, /12 and hi appear in our case. The parameter h\ can be determined by the mass difference between D 
mesons. The parameter ho is obtained by the quark mass dependence of D and Ds mesons from the lattice data [29| . 
The other constants including /13 and /15 are determined by fitting the lattice data of the scattering lengths of D and 
light mesons (24j. To confirm the existence of D* (2317), its mass and width are not used as input. All the states 
including _D* (2317) will be obtained from pole analysis on the scattering amplitudes. As a comparison, the quark 
mass dependence of Z?* (2317) is also discussed. 
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The paper is organized as follows. In sect. II, the effective chiral Lagrangian up to next-to-leading order is briefly 
introduced. We calculate the unitary scattering amplitudes and determine the low energy constants by fitting lattice 
simulations on S-w&ve scattering lengths in sect. III. In sect. IV, we present the possible bound states or resonant 
poles in appropriate channels, and then investigate the quark mass dependence of D* (2317). Finally, We make a 
brief summary in sect. V. 



II. THE EFFECTIVE LAGRANGIAN 



The leading order chiral Lagrangian for describing the interaction between the Goldstonc boson and the heavy 
pseudoscalar meson is [16l - fl8| 
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= V^DV^D^ - M D DD^ (1) 
with D = (D°, D + , D+). The covariant derivative is 

r„ = -(utfyu + ufyut) , (2) 



where 



U = exp( ^ h u 2 = Uf (3) 



with <j) containing the Goldstone boson fields 

71" T^T 



F is the Goldstone boson decay constant in the chiral limit, which we will identify with the pion decay constant, 
F = 92.4MeV. 

The strong interaction part of NLO chiral Lagrangian reads 

42. = D(-h (x+) - hix+ + h 2 (?vO - hzu^)D 

+ V^D^u^u") -h 6 {vf t ,u v }-h 6 [uf l ) u v ])VuD, (5) 

where <> stands for the trace of the 3x3 matrices, and 

X+ = U'x u + U X U ! 

Up = iu'D^Uu' . (6) 

with 

X = 2B ■ diag (m u , ma, rn s ) . (7) 

The hi term is a little different from [23j in order that the term D(x+)D will completely disappear in large Nc 
limit [19]. The corresponding coefficients C\ are also modified (see Tab. [J). The term proportional to ho leads to 
a singlet contribution to the D-meson masses which depends linearly on the light quark masses, and is the heavy 
meson analog of the pion-nucleon sigma term [30]. The hi term will contribute to the 5'C/(3)y-violating mass splitting 
amongst D mesons. The masses of D and D s mesons can be expressed as 

■ 2 

M A D = M D + Ah B{m u + m.d + m a ) + AhiBrh , 

o 2 

MS = M D + Ah B( iTiu + md + m s ) + 4:hiBm s , (8) 



from which we can determine 
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TABLE I: The coefficients in the scattering amplitudes. Here, S (I) denotes the total strangeness (isospin) of the two-meson 
system. 



where to = (m u + rrid)/2 and the mass relations of Goldstone bosons from leading order chiral expansion, 

Ml = 2Bm , Ml = B{m + m s ) , M 2 = ^B(m + 2m s ) , (10) 

are used. We can simply estimate the value of ho to be 0.055 according to the slope of the extrapolation curve from 
lattice [II]. 

III. SCATTERING AMPLITUDES AND UNITARIZATION 

The perturbative chiral amplitudes up to NLO can be easily obtained. Besides the terms in Ref. [23[ where large 
Nc suppressed ones are omitted, there are two additional terms (last two terms in Eq. Although suppressed 

in large Nc limit, the contributions from h Q , h 2 and /14 terms may not be negligible since we are working at Nq = 3. 
On the other hand, complete large Nc analysis in light meson sector shows that poles will move far away from their 
physical positions as increasing Nc |26l428|. Therefore, in this paper we include the complete tree level amplitude 
with definite strangeness and isospin up to 0(p 2 ), which can be written as 



T(s,t,u) = T w (s,t,u) +T (2) {s,t,u) 

Clo , , 2Ci 2C 35 . . 
= ^( s -") + ^^i + -^-tf35(s,t,w) 



+ ^fto + ^2r M ( a . /.,,,. (in 



where the subscripts denote the chiral dimension and the functions #35 and H24 are expressed 

H 35 (s,t,u) = h 3 p2 -pi + h 5 (pi -P2P3 -Pi + Pi -PaPi -Ps) , 
H 2 i{s, t, u) = 2h 2 p 2 ■ Pi + h 4 (pi ■ P2P3 -p4+pi ■ P4P2 ■ P3) • (12) 

Here, we adopt the same convention for the isospin decompositions as Ref. (23[. The coefficients in all the amplitudes 
are given in Tab. Q] We also dropped the h e term as in Ref. [23[ since it is suppressed by one order due to the 
commutator structure. The tree level amplitudes can be projected to the S'-wave by using 

1 f 1 

V>f' I) (s) = - d cos 6Tff' I} (s,t(s, cos 0),u(a, cos 0)) , (13) 
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where 



2 2 1 r - ' — 2 - 2ir - rni 



u(s, cosd) — m 1 + m i [s + m 1 — m 2 ][s 

2s 

-^\/A(s, mf,ml)X(s, ml, ml) cos 9 , (14) 

with 

A(s, m 2 , m 2 ) = [s - (mi + mj) 2 ][s — (m, — mj) 2 ] . (15) 

In (3lT[33j . a general method was proposed to construct scattering amplitudes satisfying unitarity, i.e. 

T(s) = V(s)[l-G(s)-V(s)}- 1 , (16) 

where V is a matrix whose elements are given by Eq. (|13[) and G is a diagonal matrix with the element being a 
two-meson integral 

n ( \ — ■ f d 4 q ii 

<Jii(, s ) — * / 77, \a ~2 2,-7 ; \2 2 , ■ ' l 17 / 

J (2ttJ* <T — m{ + te (pi + p 2 - q) - m% + ie 

with mi and mi the masses of the particles appearing in the loop. The analytic expression of Gn(s) can be expressed 
by El 



2s~ 



-{a^+log-^ + ^-log^ 

[log(s - A + a) + log(s + A + a) - log(-a + A + a) - log(-s - A + cr)]} 



(18) 



where 



o- = [-(s- (mi +m 2 ) 2 )(s- (mi -m 2 ) 2 )] 1/2 , A = ml-mj. (19) 



a(/i) is the subtraction constant with fi the regularization scale. In the numerical calculation, we tried three possible 
values of a(mr>) estimated in Ref. [2l|. Here, we did not use the the inverse amplitude method (IAM) to get the unitary 
scattering amplitude. For the light meson scattering, the divergence of one loop contributions from the leading order 
(0(p 2 )) Lagrangian can be subtracted from the next leading order (0(p 4 )) tree diagram resulting in the renormalized 
low energy constants. Since there are a lot of experimental and lattice data, the low energy constants for light mesons 
are well determined. However, for heavy meson, if we use the IAM to get the unitary scattering amplitude, the low 
energy constants at 0(p 3 ) are needed to cancel the divergence. The current data for heavy meson are not enough to 
determine all the constants at 0(p 2 ) and 0(p 3 ) very well. Therefore, in this paper, we apply the T-matrix formalism 
proposed in Refs. [3ll - l33j |. where only one parameter, the subtraction constant a(fi) was introduced. From Eq. (|2i~|) . 
one can see that G(sthr) is up to 0(p 1 ). If we include a linear dependence term of m 2 in a(/i), as pointed out in 
Ref. [lH, this higher order contribution would not change the general features of the results. 
The S'-wave scattering length is defined as 

°° = "8 7 r(Afi 1 +Af 2 ) Tthr ' (20) 

with Mi and M 2 denoting the masses of the scattered heavy and light mesons, respectively. T t hr is the unitarized 
amplitude at threshold, sthr = [M\ + M2) 2 , which can be obtained from Eq. (|T6|) . with 

1 2C 
V(s thr ) = ^ [C LO M 1 M 2 + -^hi + 2C 35 (h 3 M 2 + 2h 5 M 2 M 2 ) 

+4Coh + 4C 24 (h 2 Af 2 + h 4 M 2 M 2 )] , 

1 1 M 2 M 2 

G(s thr ) = — [a(M) + jj^M In -± + M 2 In -±)] . (21) 

There are no experimental data for the scattering of Goldstone bosons off D-mesons available. However, the low 

(2) 

energy constants entering into NLO Lagrangian £ st { can be determined by fitting the recent lattice simulations on 
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0.061 ± 0.007 


0.056 ±0.007 


0.038 ± 0.005 


h 5 
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0.032 ±0.014 


0.172 ±0.060 


Xd.o.f 


43.4/12 = 3.6 


43.6/12 = 3.6 


42.1/12 = 3.5 



TABLE II: The fit results on LECs corresponding to different values of a(mo) taken from Ref. ,21]. /12 and /13 are dimensionless, 
while /14 and /15 are in unit of GeV -2 . 
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TABLE III: The S-wave scattering lengths from calculations at LO and NLO (units are fm). The results using unitarized 
amplitudes are also given in the two columns denoted by UChPT and CUChPT, representing one-channel and coupled-channel 
unitarized chiral perturbation theory, respectively. LECs are taken from Fit III. 



S- wave scattering lengths |24j. The lattice spacing is b = 0.12fm. The s quark mass is 80MeV, which is consistent 
with its physical mass, and four ensembles are chosen with M, = 0.1842,0.2238,0.3113,0.3752 in lattice unit, or 
equivalents M w = 0.2925, 0.3554, 0.4943, 0.5958 in unit of GeV (These lattice data were misused in (23|). 
From Eq. ([5]), the pion mass dependence of D mesons up to <D(M 2 ) can be expressed as 

M D (AU) = M D \ phy + 2/ *? ± kl {Ml - MXhy) , 

MD\phy 

M Ds {M w ) = M Ds \ phy + 2h ° {Ml - MXhy) , (22) 



and from Eq. (fT0|) we can get 



Mos\phy 



M 2 o M 2 

M K {M„) = M K + , M V {M„) = M v + , (23) 

AM k 6M V 

o o 

where Mr = 486MeV and M v = 542MeV are the masses of kaon and r\ in the chiral limit, respectively. The physical 
masses for all mesons are taken from PDG [3J], i.e., MXhy = 138 MeV, Mx\phy — 496 MeV, M v \ p h y — 548 MeV, 
M D \ phy = 1867 MeV, and M D Xh y = 1968 MeV. Thus, with Eqs. ([2Tj ]> -([23" ]) . the scattering length a can be expressed 
as a function of M x . The four unknown low energy constants h,2, /13, h& and /15 can be determined by fitting the 
lattice data of scattering lengths. 

Although (S,I) = (1,1) is in fact a couple channel case, we use single channel unitarization in our fit procedure 
for all of the four channels, since so far lattice data only exist without channel coupling. The obtained LECs in 
three cases are shown in Tab. [TT1 corresponding to different values of a(m£>) which were estimated by comparing the 
dispersion relation method with the cut-off method [2l|. a{m£>) chosen to be —0.373, —0.630 and —0.864 correspond 
to the resultant cut-off momentum q max at 0.6, 0.8 and 1.0 GeV, respectively. The obtained scattering length at 
physical pion mass for each channel are listed in Tab. III. In Fig. 1, we plot the fitted results of the scattering lengths 
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FIG. 1: The fitted results on scattering lengths. The squares at physical pion mass denote chiral extrapolation results given 
by [HI, which are not included in our fit. The LECs are taken from fit III. Solid lines correspond to the central value of each 
parameter. 



versus pion mass with a(mo) = —0.864 since all the three choices of a(m£>) give similar curves. One can see that 
lattice data can be fitted very well except the (S,I) = (2,1/2) channel. The scattering lengths of (S, I) = (0,3/2) 
and (S,I) = (1, 1) channels are sensitive to pion mass which can be understood from Tab. 1. In the chiral limit, 
both of them go to zero. The scattering lengths of (S,I) = (—1,1) and (S,I) = (2,1/2) channels change little 
with the increasing pion mass. The calculation with the heavy meson chiral perturbation theory is not completely 
consistent with the lattice data of (S,I) — (2, 1/2) channel which deserves further investigation. Without the large 
Nc suppression terms, the scattering length of (S, I) = (1,1) channel remains zero with the increasing pion mass. 
In Ref. [23|], with their notation, the coefficient C\ = 2M% makes the scattering lengths of (S,I) = (1,1) channel be 
always negative and decrease with increasing pion mass. The inclusion of Nc suppression terms improves the fit of 
lattice data. Due to the inaccuracy of lattice data, we show the uncertainty of the scattering lengths in Fig. 1. The 
corresponding error bars of the low energy constants are also listed in Tab. II. For /12 and /14, the errors are about 
10% — 20% and their signs keep the same with all the three parameter sets. While for h 3 and /15, their signs could 
change with different parameter sets and the error bars are larger. 
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TABLE IV: Pole positions on y^s plane in unit of MeV. Thr and RS denote channel threshold and Riemann Sheet, respectively. 



IV. POLE ANALYSIS 



Under the present convention, the relationship between S matrix and T matrix is given by 



2z ki kj 
Sij = Sij - —^-T i:j { s ) , 



(24) 



with ki the i-th channel momentum. In general, for a system with N open channels, we have in total 2 N Riemann 
sheets, which can be enumerated as L(a%, <Ji, . . . , cr/v) where er^ stands for the sign of Imfcj [35]. Taking 2-channel case 
for example, we enumerate the 4 sheets in the following way, 



Imfci > 0, Imfc 2 > 0, L(- 

Imfei < 0, Imfc 2 > 0, L{- 

Imfci < 0, Imfc 2 < 0, L(- 

Imki > 0, Imfc 2 < 0, L(- 

The analytic continuation of 5* matrix to different sheets can be obtained 
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-iS 12 
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(25) 



(26) 



from which we can see that the poles on sheet-II and sheet-Ill correspond to zeroes of Sn(s) and detS" = S11S22 — 
S12S21, respectively. 

Corresponding to each set of parameters given by Tab. HH we list the pole positions found in appropriate channels 
in Tab. HVl from which one can see all the three parameter sets give similar results except for the (S,I) = (—1,0) 
channel. The error bars of the mass and width come from the uncertainty of the low energy constants hi. 

In (S,I) — (—1,0) channel, the pole structure is unstable. Pole positions are dependent on the strength of interac- 
tions, which is governed by LECs. In fit I, we find a "resonance", whose position is similar to [23|. However, it has no 
physical correspondence since a particle with mass below the lowest hadron-hadron threshold can not possess finite 
width by strong decay. In fit II, there is a virtual state, which is located on the real axis below DK threshold on the 
second Riemann sheet. The bound state pole predicted by fit III is in agreement with [2(|. Further experiments on 
this channel will determine which parameter set is more reasonable. 

In (S,I) = (0, 1/2) channel, we perform 3-channel unitarization. For example, for the parameter set II, we find a 
broad second sheet pole at (2122tg 1 -i93lio)MeV and a narrow third sheet pole at (2434tl-il9tV)MeV , respectively. 
Although still deviate from the experimental data 0, Q , our results are in agreement with [20, Hl| . Ref. [22[ gave 
some arguments to explain why resonances predicted theoretically have not been observed by experiment. In addition 
to production rate, finding a new state experimentally also depends on data measurements and analysis, which are 
affected by many factors, such as data statistics, the background, etc. 

A bound state pole of £>* (2317) in (5, 7) = (1, 0) channel is obtained. In Ref. [HSl], the bound state of D; o (2317) 
is assumed and its mass is used as an input to determine the LECs. Here, Z?* (2317) is really obtained from the 



analysis of the poles. For example, for the parameter set II, the mass of _D* (2317) in our analysis is m 
which is in agreement with 21 1. 



2327+^MeV, 
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In (S, I) = (1, 1) channel, only Nc suppressed terms contribute to the elastic scattering amplitudes, as can be seen 
from Tab. fl] The nearest resonance pole to the physical region is located on sheet II at (2372^25 — «39lg)MeV. The 
position of this state is different from that obtained in Ref. J21J where a resonance on sheet III with smaller mass and 
larger width exists. 

Although the mass and width of the states are obtained, we did not get the information about the nature of the 
states. There are some methods to determine the structure of a particle, as mentioned in the introduction. In Ref. [Hj], 
the authors studied the structure by investigating the quark mass dependence of the states. This meth od p rovides a 
direct and clear picture for the composition of a particle. We now make the similar analysis as in Ref. [12j . We first 
fix s quark mass and vary the light quark masses. In Fig. 2a and 2b, we show the mass of D* (2317), as well as the 
binding energy as a function of the pion mass. Both the mass and binding energy increase with the increasing pion 
mass. A pure cs state has no constituent light quarks. Its light quark mass dependence only comes from sea quark 
contributions, which should be very weak as the case of £>s(1968) shown in lattice simulations (29[. The sensitive 
dependence of light quark mass shows that Z?* (2317) may probably be a DK molecular or tetraquark state where 
the constituent light quark exists. 

Some general discussions on the importance of kaon mass dependence have been made in Ref. [12j |. The mass of 
kaon-hadron molecular state is given by 

M = M K + M h - e , (27) 

where is the mass of the other hadron and e is the binding energy. The leading kaon mass dependence of such a 
bound state is linear, and the slop is unity. 

To study the kaon mass dependence of D* (2317), we fix the pion mass at its physical value, and express all results 
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in terms of Mjc- From Eq. © and (ITUl) . we can obtain 

M D (Mk) = M D \phy + ,° (M% - M K\phy) , 

MD\phy 

M Ds (M K ) = M Ds \ phy + 2 ^ + , kl) {M 2 K - M 2 K \ phy ) , (28) 

JMDs\phy 

and M V {M K ) = ^\M\ - \Ml\ phy . 

Fig. 2c and 2d show the mass and binding energy of Z?* (2317) as a function of kaon mass. The if-meson mass 
dependence of D* (2317) is almost linear which is really in good agreement with the DK molecular expectation. 
These results are comparable with Refs. (l2l |36||. 



V. SUMMARY 



In this paper, we calculate the complete scattering amplitudes of Goldstone bosons off the pseudoscalar D-mesons 
using unitarized heavy meson chiral approach. Two low energy constants ho and hi are determined by the mass 
splitting among D mesons and pion mass dependence of D and D s . The other four LECs are determined by fitting 
lattice simulations on S'-wave scattering lengths. The large Nc suppressed terms improve the fit. Three sets of 
parameters are obtained according to different choices of the subtraction constant a(mo). All the three sets of 
parameters give similar scattering lengths which are close to the lattice results. 

For three parameter sets, the positions of the poles in each channel are close except in (S,I) = (—1,0) channel. In 
this channel, the poles are sensitive to the parameter sets. Further experiments or lattice simulation can determine 
which parameter set is more reasonable. For other channels, we take the average values of the mass and width as the 
final results. Z?* (2317) is obtained as a bound state in (S, I) — (1, 0) channel, with the mass being m = 2326t^MeV. 
The strong pion mass dependence of its mass and binding energy disfavors conventional cs content. The approximately 
linear kaon mass dependence reveals it is predominately a DK molecular state. In (S,I) — (0,1/2) channel, a 
broad pole structure is found at (2124+Jq — i97*\\)MeV on the second Riemann sheet, and a narrow pole is at 
(2429±io - H5±^)MeV on the third Riemann sheet. A resonance pole also exists on the second Riemann sheet with 
mass and half width 2374±g MeV and 34±^ MeV in (S,I) = (1, 1) channel. 
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